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IMPACT OP A VEE-TYPE SEAPLANE ON WATER WITH 



REFERENCE TO ELASTICITY* 
By F. Weinig 



The theory developed "by H. Wagner for the computation 
of the landing impact on water for a rigid float is extend- 
ed to include elastic floats "by introducing the concept of 
an equivalent rigid bottom to substitute for the actual 
elastic bottom. 



OBJECT OF INVESTIGATION 



The theory developed by H. Wagner for the computation 
of the landing of float bottoms on the water is extended 
to include bottoms having elasticity in order to take ac- 
count of the elasticity factor on the landing impact. 



WAGNER'S THEORY OF THE IMPACT OF VEE-TYPE FLOATS 



WITH RIGID BOTTOMS 



Let the downward velocity of the seaplane V just 
before impact be denoted by V Q . During the immersion 
process which is assumed to start at t = 0, V = V(t). 

In investigating the impact we assume for convenience 
that the float is at rest while the water is in motion 
relative to it. Essentially this flow is similar to that 
of an infinite fluid about a flat plate at rest, the width 
2c of the plate corresponding to the instantaneous width 
of the impact area of the float (fig. !)• T ke velocity of 
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the water particles at the free surface (2 > c) is 

V 



v = v 



n - - , g- 



A- 



c 



The rise Y of the water measured from the instant of im- 
mersion is 

t t 7 dt 

Y = / v n dt = / 



'J 'B. - J I 

o o / i ^ £ 

v X 



5" 
X 2 " 



The width of the impact area increases with time; i.e., 
c = c(t), so that since t = t(c), V = 7(c) and with 



c~x 

, \ 7 5 u(c) dc 

, a -° J 1 - 

As soon as the water particle, at position X reaches the 
edge of the impact area, c = X and Y = Y-^ , sc that 



*li - J == 



We set 

y = Y * • x - x • t - c ■ ti - r l u(n H (i) 
y 172' " 172' 1 - 172' 71 V / , t \* 

J 1 " (I) 

Let the "bottom shape be expressed ty : 

Tl(x) = p x x + p 2 x a + 0 3 x 3 + '. . .-. (2) 

Then u( \ ) = Yi + Y 3 | + Y 3 | 3 + 

c £ 

Since c < X, if we set — = = sin a, then expression 

A X 

(2) becomes T] ( x ) = Pi i + p 2 i 2 + P3 x 3 .... 
■rr/2 

= x / (Yi + Y 2 x sin a + Y 3 x 2 sin 2 a + ....) da, 
0 
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u( |) = kj. px + k 2 Pa 6 + k 3 I + .... 



(2a) 



where 
ki = 

k 3 = 

k R = 



2 2 
- X — 

TT 1 



£ X £ x i 

TT 1 3 



=0.636; k a = 1 



= 1.272; k 4 = § 

= 1.696; k 6 = § X £ 



1.000 
1.500 
1.875 



^(3) 



= £x2_ x H x £ = 2.0U0; k s = I x ^ x ^ = 2.18 



In general , 



k n = 



tt/2 

J sin 11-1 a da, 



From k n there is derived the formula 



k n +i = 



2n 



TT k 



n 



For large values of n the following approximate for- 
mula applies : 



while for small values of n the approximate formula is 

(n < k) 



2n - H + 



TT 



TT 



In some cases the "bottom shape is "better expressed "by 
T| = p x x + P n x n 
and the value of u is then obtained as 

u = k x + k n p n x 11 " 1 
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For nonintegral values of 
"be read off figure 2. 



n the values of k r 



may 



In order to determine the force P exerted on the 

float, we again consider the surface of the water to "be at 

rest and the float moving relatively to it with velocity 

V = V(t). The momentum of the fluid at the lower half 
olane is known to "be 



J=|pc 2 V=n.p 



dc V 

from which, substituting jr = — , we obtain 

du u 



p = 



dJ 
dt 



- TT 



" (I; 



u 



dV 
dt 



(5) 



The momentum of the seaplane of mass m i s equal to 
m (V 0 - V) and for the case we are here considering where 
the elasticity, for example, of the landing gear and float 
bottom is neglected, this is equal to the momentum impart- 
ed to the fluid. With 



U = 



2m 



(6) 



and 



5 p (!) V g = m (V Q - V) 



V = 



i + v. V 



(7) 
(8) 



Prom expression (S) it may be seen that V( | ) does not 
depend on the bottom shape-but only on the width of the 
wetted surface and on the mass. If the work of deforma- 
tion is neglected then 

dV 



P = - m 



dt 



(9) 



and we thus obtain 



P = TI p 

Por briefness, we set 



2 ° (1 + ^ 2 ) 3 u(| ) 



(10) 
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so that 



Pi 



= B 



B = 



(i + m. r) 



(id 

(12) 



The pressure di stribution on the float bottom will now 

2 2 

he determined. For X < c the velocity potential is: 



cp = - v/c 2 - x 2 = - v c y\ - ( 



The fluid pressure p is, in general, 
P = - P 



In the case ^ e are considering F(t) is constant 
since there is no superposition of external variable pres- 
sures. On the surface of the water F(t) = - p^/p. If p 
is the pressure above the air pressure p^ then 



P = - P 



We have 



+ i v 2 

3t 2 x 



^x = - 
x 3X 



so that 



E = L 
p u 



2 



dt J W 2 



2 - 1 



(13) 



The last term of the above equation does not apply to the 
edge of the impact area and is small compared to the first 
t\?o terms. 



P =. 



With the aid of (8), (9), and (10), we obtain: 

2 '' " u(l + M- & 3 ) 



.ft 



LULL 



-2n i' 



-(f)- 



(I) 2 - 



(1U) 



This formula does not give correct values for the edge of 
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the impact area. There the flow is similar to that about 
a planing surface moving sideways with velocity 



1 dt u 

The maximum pressure that is set up is therefore, 



P max = *i = I V i 2 



(15) 



(16) 



or with 



. 1o 



max 



,.3 s ^ .2 



(17) 



THE CONCEPT OP AN EQUIVALENT RIGID BOTTOM TO 
REPLACE THE ELASTIC BOTTOM 



Prom equation (12) it may "be seen that B_( | ) depends 
only on the form parameter of the float bottom. Such a 
parameter may also be found for a float with elastic bot- 
tom. Let y re i(|) De the elastic displacement at position 
| due to the load P(£)» then 



n(|) = T1 0 U) + T| re i(£) 



(is) 



where T] 0 (|) represents the shape of the nondeformed 
float bottom and T| ( % ) the shape of a nonelastic float 
bottom which replaces that of the elastic float. The re- 
sults obtained for the rigid float may thus be applied to 
the case of the elastic float. The equivalent rigid float 
fully takes into account the rate of expansion of the width 
dc/dt due to the elasticity; that is, the main part of the 
deformation; and only a small part of the deformation is 
neglected, namely, that of the bottom opposite the chord 
between the keel and edge of the impact area. 

We must next consider how y re i is to be determined. 
It may be seen from equation (lH) for the pressure distri- 
bution that the relative distribution of the pressure is 
essentially independent of the shape of the float bottom, 
i.e., of u( | ) since the last term in the brackets is very 
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small compared to the other terms for the greatest part of 
the range 0 < x < |. 

If the elastic deformation for the same relative pres- 
sure distribution is proportional to the total load, then 
the deformation under the load ■ B ( \) at any position | 
may be expressed in the form 

Tl r el(£) = 1(1) (a 0 + ai £ + a 2 f + a 3 f + ) (19) 

or with I = x. 

Tl rel (x) = B_(x) (a 0 + a 1 x + a 2 x 3 + a 3 x 3 + ....) 



in parenthe- 
A(e): 



For briefness we shall denote the expression 
ses which takes account of the elasticity by 

(a 0 + ai | + a 2 | 2 + a 3 | 3 + ....)= A(£) 

The coefficients a 0 , a x , a 2 , a 3 must be determined, 
for each type of float, either by experiment or computa- 
tion. The pressure distributions given by equation ( lU ) , 
which applies to rigid floats, are to be substituted and 
thus the bending at the edge of the impact area determined'. 



We have 
port i onal t o 
of the float 



assumed that the elastic deformation was pro- 
the load P, but due to the small thickness 
bottom this assumption will in many cases not 
prove to be correct. In these cases, however, it may be 
taken as a first approximation that for a moderate in- 
crease in an assumed loading the increase in the deforma- 
tion is proportional to the load increase. 



In order to determine the deformation it may be as- 
sumed first that the load is the same as that exerted on 
the nondeformed float on impact, the corresponding defor- 
mation under this load determined, and. then the deforma- 
tion when the load is slightly increased. In this way an 
expression for the deformation may be found of the form 

Tlrel( x ) = lorel^^ + a o + a x x + a 2 x H + ....) 
If we then add T| 0 and 11o r el' we °' b ' tain: 

"Ho + Tlorei = ^o = Pi s + ^2 x 2 + p 3 x 3 + .... 
And the same formulas apply for the determination of the 



s 
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impact force as for the case where the deformation was 
taken proportional to the load, provided we replace Pi 
by Pi. 

This method of computation is admissible provided the 
equation for y re i gives sufficiently accurate values 

corresponding to the load P first obtained; otherwise, 
the computation must he repeated for a new value of P. 

PROCESS FOR OBTAINING THE EQUIVALENT FLOAT BOTTOM AND 
DETERMINING- THE EFFECT OF THE ELASTICITY 
ON THE IMPACT FORCE 



a) Mathematical f ormulation of the problem as an in- 
tegral equation .- We- have given Tl 0 ( £ ) , p., A(|) and seek 

to determine 1"l re l» From equations (12) and (ig) 

Tlrel = 5 ~ (20 > 

(1+m. I s ) (u 0 (|) + ureld)) 

and therefore, 

u(|) = u 0 (|) + u rel (|) = LMAl 

(l + n t 2 ) Tlrel 

and from equations (l) and (IS) 

so that for the determination of "n re l' we o* 5 * 3 -* 11 tntJ 
following relation: 

This is a nonhomogeneous , nonlinear integral (Volterra) 
equation of the second class. No direct methods have so 
far been developed for the solution of equations of this 
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kind. In the case, however, where there exists a solution 

u(|) = 0 the iteration method may probably lead to the 
solution. 

"b ) Development into series and approximate solutions . - 
Instead of carrying out the iteration process a solution 
was first attempted "by development into series and the 
following series were'tried: 

A. " = P Q + p x | + p 2 | 2 + . . . 

B. ir (1 + n f f = p o i + Pi > i + p 2 ' | 2 + . . . 

°' ^rel = Pl ' x + Pa ' xS + p s 1 x 3 + ... 

These expressions may he simplified, depending on the type 
of float, hy the elimination, for example, of the first 
term of the series. None of these series, however, con- 
verges in general (corresponding to the non con verging of 

the series in the development of = l- x + x 2 -x 3 + ~ 

... for x > 1 ) . 

It was then attempted to find approximate solutions 
hut it was not possible to ohtain good approximations since 
the system of equations to he solved was no longer linear 
so that this method, too, was unsatisfactory. 

The possibility of using the iteration method was not 
tried since it could only he applied to individual cases 
in connection with a graphical process, and the results 
thus obtained would then have no general significance. 

c ) Obtaining solutions by applying the inverse proc- 
ess .- To obtain results of more general application, the 
inverse process may be employed; that is , we start out 
with the deformed bottom shape and seek to find the origi- 
nal undeformed shape. 

We thus have as given T] = T\ Q + T| rel and the function 
A(|), and we obtain in turn P, T| re i. and T| 0 : 

ti 0 = nee > — a(|) (22) 

(i + n e 2 ) u(i) 



As an example, we choose the values 



■ Ill IB mm* I !■ ■■INMHHIIII II 



■ ■II III ■IIIIBBHHM ■■IIMIM I I II I I III I I III I I !■ I II III II I I 
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P, = 3 X x 

M. = 0 

A(|) = cp R(|) = (| P X ) 2 * R(|) 

R ( t ) = a 0 + a x | + a 3 £ 2 .... 



(23) 



which represent a straight V "bottom float of a seaplane 
having a very large mass. We thus obtain: 



~ cp H(x) 
- 3i 

TT 



or with 
cp 



Yi' 



Tl 0 = Tl (l - | \|/ R(x)) = p x (l - | M/ a 0 ^) x 

- p £ \|/ (a,! x 2 + a a x3 + a_ x4 + .,.) 

TT 

ki P x ( 1 - £ \J/ - p x 2 \|/ (k 2 a x x + k 3 a 2 x 2 + ...) 

n TT / TT 



Un = 



B = 



2o = 



We have 



kx Pi 



kiPi (l - - i a 0 j -pj_ £ v|/(k s a! |[ +k 3 a 2 | 2 + ) 

\ TT / TT 



so that 



B 
% 



= 1 - v|/ (k x a 0 + k 2 a x % + k 3 a 3 | 3 ) 



and at | = 1, where the maximum impact force is to be ex- 
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pect ed , 

Irel _ 1(1) 



— O — c 



j^-j - 1 = - >j/ (ki a 0 + k 2 tti + k 3 a 3 



The deformation is: 

Tlrel = B(|) A(f ) = q> H ( | ) 

and attains its maximum value for 

<£__ [R(n + | R'(£)] = 0 
ki Pi — — 



From the above equation 

R 

i'U) 



I = I = - ^-4- (if I < 1) 



and therefore 

cp R 2 (f) 



r el 



or 



T) (1 ) 

_ ^rel 1 R 1 ( | ) 

" TKiT k7 R 2 (!) 

and thus ( for | < 1 ) 

For | > 1 the greatest deformation is at 
so that 

n rpl (D 

v ~ Ticrr ki r(d 
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and therefore, 

jjrel _ _ 'HrelC 1 ) tt k x a 0 + k a + k 3 otg + ... 
B_ 0 T] ( 1 ) 2 a 0 + a x + a 2 + . . 

As a first example (fig. 3a), let the following val- 
ues "be chosen: 

11 = 01 x 
IX = 0 

and A = cp 1(1-1) 

so that 

R = |(1 - i) (a 0 = 0, a x = 1, a 3 = - 1) 

In this case the bottom (construction A) is deformed be- 
tween the keel and the chine in such a manner that for a 
constant impact force Tlrel would be a parabola. 

The maximum deformation is at 

I=-_(1_Z_I) ; T= 2 
1 - 2 | 3 

so that 

~£el = £7 (U - „) W = 2.S 9 5^ 
B Q g 71(1) Tl(l) 

The relative increase in the impact force in this case 
is therefore 2.89 times the ratio of the maximum deforma- 
tion 'Hrei to the dead rise T|(l). 

As a second example (fig. 3b), let the following val- 
ues be chosen: 

Tl = $i x 

p. = 0 

and A = cp £ 2 

so that 

R = 1° (a 0 = 0, a x = 0 , a 2 = 1 ) 
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In this case (construction 3 "b ) , the bottom is elastic 
at the sides so that the chine moves up on impact. The 
form of the .A function is then obtained by assuming a 
rigid bottom with the impact force approximately concen- 
trated near the edge of the impact area. 

The maximum deformation is now at x = 1, so that 
Irel _ _ Tl r ei(l) ' 



B 0 Tl(l) 

This construction therefore leads to a decrease in 
the impact force, namely, by twice the ratio of the defor- 
mation at the edge to the final dead rise T| ( 1 ) . 

For the third example (fig. 3c), let the values chosen 

be : 

n = Pi x 
n = c 

and A = cp ( 1 - | ) 



so that 



R = 1 - | 2 (a 0 =1, a x = 0, a s = - 1) 



(construction C). 

The maximum deformation is at 

I = 1 - F. T - /I 



2 i 



3 



so that 



Irel _ 3 /s-^rel _ ? c ^rel 

IT " a vB TUTT " 2 * b TfUT 

In this case, therefore, the relative increase in the 
impact force is 2.6 times the ratio of the maximum defor- 
mation to the dead rise. 

From these examples it may be deduced that under cer- 
tain conditions the effect of the elasticity of the bottom 
may be of some importance but in especially unfavorable 
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cases the effect may even be more considerable, so that no 
conclusions of general application can "be attached to the 
results obtained from the examples just given. 



SUMMARY 



For several particular cases of possible float-bottom 
shapes the impact force due to the elasticity deformation 
was compared with that obtained for the rigid bottom. The 
work carried out considers, however, only part of the ef- 
fect of the elasticity on landing and obviously the less 
important part. A much more important effect is that of 
the elastic connection between the principal masses of the 
seaplane and the floats, and this problem still awaits so- 
lution. 



Translation by S. Reiss, 
National Advisory Committee 
for Aeronautics. 
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Pigs. 1,2,3 



y'f n 




X = B/2 
c = B/2 

Figure 1.- Notation for rigid V-type float bottom landing on 
water . 



2.0 



1.6 



1.2 



.8 
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(a) HI 





n 



Figure 2.- Values of k n> 



II : 

(c) Ilir 
Figure 3.- Several types of 

possible deforma- 
tion for float bottoms, 
(a) Points 1,11,111,1V, and V 
are mutually fixed and 
bottom is deformed between II-III and III-IV. 

(b) Points I, IIJ and V are mutually fixed and bottom is deformed at 

II and IV. (c) Points I, II, IV and V are mutually fixed and deformation 

is possible at the keel III. 
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